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ABSTRACT 

We discuss interpretation of the cosmic infrared background (CIB) 
anisotropies detected by us recently in the Spitzer IRAC based measurements. 
The fluctuations are approximately isotropic on the sky consistent with their 
cosmological origin. They remain after removal of fairly faint intervening sources 
and must arise from a population which has a strong CIB clustering component 
with only a small shot-noise level. We discuss the constraints the data place on 
the luminosities, epochs and mass-to-light ratios of the indvidual sources produc- 
ing them. Assuming the concordance ACDM cosmology the measurements imply 
that the luminous sources producing them lie at cosmic times < 1 Gyr and were 
individually much brighter per unit mass than the present stellar populations. 

Subject headings: cosmology: theory - cosmology: observations - diffuse radiation 
- early Universe 

1. Introduction 

If the early Universe contained significantly more luminous populations than at present, 
such as is thought to be the case with the very first metal-free stars (see review by Bromm 
& Larson 2005), these populations could have produced a significant contribution to the 
cosmic infrared background (CIB) with potentially measurable structure (Santos et al 2002, 
Salvaterra & Ferrara 2003, Cooray et al 2004, Kashlinsky et al 2004; see Kashlinsky (2005) 
for recent review). In an attempt to uncover the CIB fluctuations from early populations 
we have analyzed deep images obtained with the Spitzer IRAC instrument (Kashlinsky, 
Arendt, Mather & Moseley 2005; hereafter KAMM1), which led to detecting significant CIB 
fluctuations remaining after subtracting sources to faint flux levels. In a companion paper we 
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presented analysis from deeper and larger fields using the GOODS Spitzer data (Kashlinsky, 
Arendt, Mather & Moseley 2006; hereafter KAMM2), which confirms our earlier findings 
and extends them to fainter levels of removed galaxy populations and larger angular scales. 

In this Letter we discuss the cosmological implications of the recent measurements of 
the CIB fluctuations from early populations obtained by us (KAMM1,KAMM2). These 
measurements imply that the signal must come from cosmic sources which have a significant 
clustering component, but a low shot-noise contribution to the power spectrum. Given 
the amplitude of the CIB flux expected from these populations in the concordance ACDM 
cosmology (>1 nW/m 2 /sr ), we show that these sources must have very faint individual 
fluxes of < 10-20 nJy in order not to exceed the measured levels of the remaining shot-noise. 
Furthermore, these populations had to have mass-to-light ratio significantly below that of the 
present day stellar populations in order to produce the required CIB fluxes in the short cosmic 
time available (<1 Gyr) from the available baryons. Finally, we discuss the prospects for 
their individual detection with future space missions. We use the AB magnitude system, so 
flux per frequency v of magnitude m is S u (m)=3631 x io~°- 4m Jy; diffuse flux in nW/m 2 /sr is 
defined as vl u , with I u being the surface brightness in MJy/sr. 



2. Magnitudes and epochs of the sources of the CIB fluctuations 

In their analysis KAMM1 and KAMM2 used a total of 5 different fields with deep 
Spitzer IRAC observations of up to 24 hours per pixel. All the observed fields are located at 
high Galactic and Ecliptic latitudes and are free of significant zodiacal emissions at all IRAC 
channels and of cirrus at the IRAC channels 1-3 (3.6 to 5.8 /iin). Individual galaxies and 
other sources were removed until a fixed level of the shot noise from remaining sources was 
reached. The power spectrum of the remaining diffuse emission showed a residual shot-noise 
component on small angles and a significant excess due to clustering of faint / distant sources 
at scales ^0.5'. Within the errors all fields cleaned to the same shot-noise level showed the 
same excess fluctuations consistent with their cosmological origin (see Fig. 1 of KAMM2). 
At 8 /iin there is a significant pollution by the Galactic cirrus and at 5.8 fim the larger 
instrumental noise leads to relatively large errors in the large-scale fluctuations. Here we 
concentrate on the interpretation of the data at 3.6 and 4.5 /xm in terms of the luminosities, 
the epochs and the nature of the cosmological sources contributing to these fluctuations. 

KAMM1-2 show that the CIB fluctuations must come from cosmological sources, such as 
ordinary galaxies and the putat ive Population I II. The former are defined as metal-rich stars 



with IMFs of a Salpeter-Scalo (IKennicutt 19981 ) type with masses ~ 1M . Population III is 



defined (loosely) as luminous sources that existed at, say, z> 10 which possibly were individu- 
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ally very massive and intrinsically very luminous. Data such as discussed here cannot resolve 



whether the sources contributing to the CIB were metal-rich ([Fernandez Komatsu 20051) 



and whether the source of this radiation was stellar nucleosynthes is ( ISantos et al 20021 ) or 



black-hole accretion in the early Universe (jCooray &: Yoshida 2004). Population HI epochs 



z > 10, may contain emissions by both stars and quasar-like objects ( iKashlinsky fc Rees 19831 ). 



Any model aimed to explain the CIB fluctuations results must reproduce three major 
aspects: 1) The sources producing the measured CIB fluctuations must be fainter than those 
removed from the data. 2) They must reproduce the observed excess CIB fluctuations at 
>0.5', where 5F ~0. 07-0.1 nW/m 2 /sr . 3) Lastly, the populations below the above cutoff 
must account not only for the correlated part of the CIB, but must also reproduce the 
observed (low) shot-noise component of the signal. These lead to: 

1) The the shot-noise component of the power spectrum from source counts dN/dm 
per magnitude interval dm is Psn= / S 2 (m)dN(m) with S=vS v . To estimate limiting mag- 
nitudes implied by the measured shot-noise, w e generated source counts for the observed 



fields with SExtractor (IBertin fc Arnouts 19961 ). Fig. [T] shows the remaining shot noise 



levels in KAMM1-2 analysis and the counts data. The intersection of the counts with the 
lowest shot noise levels shows that the sources are eliminated to mt 25-26, so the detected 
CIB fluctuations come from fainter sources. This magnitude limit at 3.6 /im corresponds 
to only lO 9 /i~ 2 lO~ a4(m - 25 - 5) L emitted at 6000 A at z=5 where h is the Hubble constant 
in uni ts of 100 km/sec /Mpc. If the counts contain extra populations in addition to those 



from flFazio et al 20041 ). the magnitude limit will be fainter. Thus KAMM1,KAMM2 have 
removed a significant fraction of galaxies even at z=5 and the CIB fluctuations must come 
from sources at higher z. 

2) The clusterin g compon ent of the CIB at 0.5'^27r/g<5 / requires Fqib ~ afewnW/m 2 /sr as 



noted by us earlier (IKAMMll ). The rms fluctuation in the CIB flux, 5F=^Jq 2 P 2 (q)/2'K, on 
angular scale 2ir/q is related to the underlying 3-dimensional power spectrum of the emitters' 
clustering, Pz(k), the duration over which the flux was produced, At, and the rate of the 
CIB production rate, dF/dt, via the Limber equation (e.g. Kashlinsky 2005a): 

t\F F A -. 2 _ ^UdF/dt)'AHqd-/)dt 



where A(k) = {k 2 P 3 (k) ^ircAt] 1 / 2 is the rms fluctuation in source counts over the cylinder of 
radius 2ir/k and length cAt. In the limit when the CIB release rate is approximately constant, 
the relative CIB fluctuation, A F , will be ~ (A 2 (qd A 1 )) 1/2 with (...) = (At)- 1 J Af ... dt. If 
dF/dt peaks at some cosmic epoch z p , the relative fluctuation will be ~ A(qd~J^ (z p )) . 



To evaluate the range of the expected CIB flux from the sources producing the mea- 
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sured fluctuations, we adopt the ACDM model with (O, f2b ary on, ^a, /i)=(0. 3, 0.044, 0.7,0. 71) 
and consider the epochs spanning 5< z <20. The cosmic time at z=20 is ~0.2 Gyr and 
the time between 2=20 and z=5 is 1 Gyr. The scale r 8 =8/i _1 Mpc, with today's density 
contrast a 8 , subtends 6*8 ~(3-4)'. The relative fluctuation in the projected 2-dimensional 
power spectrum, A, on that angular scale 9$, produced from sources located at mean 
value of z and spanning the cosmic time At, would be A(0g) ~ <7g(l + z)" 1 ^/ cAt) 1 ' 2 ~ 
0.02a 8 (^)) _1 (^:) _1//2 neglecting the amplification due to b iasing. Bias ing, due to sources 
forming out of rare peaks of the density field, will increase A (jKaiser 19841 1 and for reasonable 
bias factors (~ 2 for systems collapsing at z ~5 to ^3 at 2^10) one can gain amplification 
factors, A in A o f ~ 2 to 1 4-5 between z=5 and 20 (IKashlinsky 199ll ; iKashlinsky 1998 



Cooray et al 20041 ; IKashlinsky et al 20041 ). Thus the arcminute scale CIB fluctuations of 



5F ~0. 07-0.1 nW/m 2 /sr at 3.6, 4.5 /im require the mean CIB from these sources to be 
F C m ~ 4 - 5(A[(1 + ^)/6]" 1 )" 1 (At/lGyr) 1 / 2 nW/m 2 /sr . Assuming that the fluctuations 
are produced by low surface brightness systems at much lower z does not alter the required 
high value of their mean CIB contribution: e.g. taking At=5Gyr corresponding to the 
cosmic time between z=l and 20 gives A ~0.02 at V assuming no biasing. (As discussed 
below such sources would likely produce shot-noise in excess of what we measure.) We 
can reach similar conclusions with the entire range of scales ^0.5' where we measure the 
clustering component of the CIB. The left panels of Fig. [1] show the least squares fits for 
Fcib, assuming the ACDM model, from all the fields data at 3.6 and 4.5 /mi. This gives 
F CIB (A(i T t£))>(4,2.5)(At/lGyr)- 1 / 2 nW/m 2 /sr at (3.6,4.5) /an respectively. 

We thus conservatively take the fiducial flux of -Fcib=1 nW/m 2 /sr as the minimal CIB 
flux at 3.6 and 4.5 /mi required by the fluctuations, corresponding to the relative mini- 
mal CIB fluctuations of ~7%. The results below can be re-scaled to arbitrary Fqib, but 
our general conclusions will be valid unless the CIB flux from sources producing the mea- 
sured fluctuations is significantly below the above number. Although the net CIB fluxes 
may be, in principle, much higher, this minimal CIB level at 3.6 /im is smaller than the 
claimed CIB excess from DIRBE and IRTS measurements over th at from galaxy counts 
JPwek fc Arendt 19981 : krendt fc Dwek 20031 : iMatsumoto et al 2005h . and is consistent with 
the recent measur ements of absorp t ion in the spectra of fa irly distant (z=0.13-0.18) blazars 
at TeV energies (IDwek et al 20051 : lAharonian et al 20061 ). Such CIB levels should, how- 
ever, be measurable from the spectra of gamma-ray bursts at z>l-2 de tectable with the 
upcoming NASA 's GLAST mission out to 300 Gev (IKashlinsky 2005bl ). Spitzer counts 
( jFazio et al 20041 ) show that the remaining ordinary galaxies can contribute only ~0.15 
nW/m 2 /sr at 3.6 /mi ( IKAMMll ). so to explain the CIB fluctuations with the remaining 
(extrapolated) Spitzer counts sources requires almost ~100% fluctuation on arcminute scales. 



3) The CIB in the populations producing the measured fluctuations significantly ex- 
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ceeds that from extrapolated IRAC counts (IFazio et al 20041 ). so the excess flux must come 



from faint er populations with a significant deviation from the extrapolated counts slope 



( IKAMMll ). The measured fluctuations indicate a population with a relatively strong cluster- 
ing component, which at the same time has low shot noise. This means that the sources must 
be individually faint. The shot-noise from the remaining galaxies dominates the power spec- 
trum of the CIB at <0.5' and its amplitude sets an upper limit on the shot-noise component 
of the sources contributi ng to the arcminu te scale CIB fluctuations. The amplitude of the 



shot-noise component is (IKashlinsky 2005al ): Psn= J >m S(m)dF(m) = S(rh)F tot (>m), where 
dF(m)=S(m)dN(m) is the CIB from sources at the magnitude interval dm and F to i{m) is 
the total flux from the remaining sources of > m. The sources contributing to the clustering 
component of the fluctuations at arcminute scales must not exceed the level of the residual 
shot noise in the data of Psn — (2, 1) X 10~ n nW 2 /m 4 /sr at (3.6, 4.5) /im. At 4.5 fim this 
shot-noise amplitude of PgN=10~ n nW 2 /m 4 /sr or 10(A/3/zm) _1 nJy-nW/m 2 /sr , would lead 
to sources contributing to the signal having mean fluxes less than 12 (FciB/ n Wm~ 2 sr _1 ) _1 
nJy or AB magnitudes fh > 29 + 2 . 5 lg ( nWl ^ -^r - 1 ) • At 3-6 /im the shot-noise levels are 
a factor of ~ 2 larger leading to m about one magnitude brighter. An important further 
information could be obtained in still deeper measurements by setting a lower limit on the 
shot-noise component of the sources contributing to the CIB fluctuations determined when 
the clustering component disappears of is substantially reduced. 



3. Discussion 

More information on the nature of populations of these faint sources can be obtained by 
considering the fraction of baryons that went through stars prior to z 1 5 (At < 1 Gyr) needed 
to explain the level of the CIB required by our data. The net flux at frequency v produced 
by the population with comoving luminosity density £, is Pcib=^ f A C U >(1 + z)~ 1 dt, where 
v'=v{\ + z). This requires the average comoving luminosity density at (0.36-0.45)/zmy^ of: 

l * ^F CIB (Atr(l + ^1.2x lO°L Mpc-s ffil±£ F ™ ; (2) 
c At 10 nWm sr _1 

For comparison the present-day luminosity density measured by t he Sloan Digital Sky Survey 



at 0.32 /im to 0.68 /im is about an order of magnitude lower (IBlanton et al 20031 ). This 
indicates significantly more luminous populations contributing to the CIB fluctuations than 
at present. The contribution to the density parameter by these sources is thus given by: 

( r £) I (0.36-0.45)^^ _ 3 F C m r / At V 1 1+2 

* = 3H$/8nG L ^ 8 - 3x10 nWm-^sr-i ^ ) ~W (3) 
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where T is the mass-to-light ratio. Fo r comparison the mean density in pr esent day stars is 
significantly lower at f2* jZ=0 ~ 2xl0~ 3 ( IFukugita et al 19981 ; ICole et al 20011 ) and much of the 
contribution to fi* iZ= o comes from the late stellar Population I stars with solar metallicities. 
Strictly speaking eq. [3] assumes no re-processing of baryons and may thus overestimate 
the required amount of luminous baryons in the case of short lived massive stars, such as 
Population III, but it shows that it is energetically easier to produce the significant CIB levels 
implied by the Spitzer data in the cosmic time available with stars whose mass-function is 
skewed toward T -C r . (For populations made up of massive stars it can be replaced with 
eq. 3 of Kashlinsky, 2005b). If the CIB fluctuations are produced by populations containing 
a significant fraction of low-mass stars, which should still be burning today, they would 
require a large fraction of the present-day stars to have been produced at z>6-10. 

To model the ordi nary stellar populations, we hav e run stellar evolution models us- 
ing the PEGASE code ( jFioc fc Rocca-Volmerange 19981 ) . assuming normal IMF with vari- 
ous metallicities and the ongoing star formation (i.e. star formation rate oc exp(— t/tburst) 
with thnrst=20 Gyr ). For Population III we adopted the spectral energy distribution from 
( iSantos et al 20021 ). Fig. [2] shows the luminosity per unit mass in stars (r -1 ) assuming the 
ordinary population to be less than 1 Gyr old (r ~O.2-O.5r ) and contrasts them with the 
expectations for massive Pop III systems (lO _2 -lO _3 r ). If the CIB fluctuation signal comes 
entirely from the Population III systems, eq. [3] would give the minimal fraction of baryons 
locked in them ~ 0.15%(Fcib/iiW m _2 sr _1 ). If the baryons are re- used in stars this fr action 
would be decreased. This number is in agreement with that of (IKashlinsky 2005bl ) after 
scaling to the appropriate CIB levels: 0.14% (F CIB ,boiomctric/nW m" 2 sr~ 1 )(z/10)(e/0.007) _1 
assuming the hydrogen burning efficiency e. (Such massive stars would be fully convective 
with the overall efficiency reaching e>3 x 10~ 3 , Schaerer 2002). 

The sources satisfying the above constraints had masses in luminous matter of: 



M* ~ 4tt4(1 + z)~ l TS(m) < 7 x 1O 5 /T 2 M 



(3.6-4.5)Mm 



[S(m)/v] fl + z 



5 x lO" 3 r 20nJy 



10 



1.6 



(4) 



where the luminosity distance was approximated di — 3.2(1 + z) 1,3 h ^pc. Such Pop 
ulation III sources, with only < a few times 1O 5 M in stellar mat erial, would be be 



low the detection thres hold in the high -z Lyman dropout searches of (IBouwens et al 2005 



Willis fc Courbin 20051 ) considered by (ISalvaterra fc Ferrara 20061 ). In any case theoreti- 
cal predictions of the luminosity function of Population III sources are necessarily model- 
dependent as they depend on the assumptions of the small-scale power and its evolution as 
well as the microphysics governing the various feedback effects during the collapse of the first 
haloes. The Press-Schechter type presc riptions may break down for the slope and regime 
of power spectra on the relevant scale (ISpringel et al 20051 ) and the feedback mechanisms 
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due to the H2 destruction by the Lyman- Werner bands radiation (jHaiman et al 19971 ) likely 
suppress star formation in a complicated halo-mass dependent way. 

To resolve the faint sources responsible for the CIB fluctuations, their individual flux 
must exceed the con fusion limit usu ally taken to be a > 5 times the flux dispersion produced 



by these emissions ( iCondon 19741 ) . Lower flux sources will be drowned in the confusion 
noise; of course, this is precisely where CIB studies would take off. If such sources were to 
contribute the CIB required by our data, at 3.6 /im they had to have the average surface 

density of n ~ Fcib/ p sn ~ 2 arcsec" 2 ( nW ^ B sr -i ) ( 10-n nwffn-^-O • To avoid the 
confusion limit and resolve these sources individually at, say, 5-sigma level (a=5) one would 
need a beam area Wbeam < a~ 2 /n ~ 0.017 ( ~uF^pr ) arcsec 2 or circular radius <0.07 
(FcrB/nWrn~ 2 sr _1 )~ 1 arcsec. This is not in th e realm of the curre nt instruments, but the 



JWST could be able to resolve these objects (IGardner et al 20061 ). Extrapolation of this 
argument to shorter A is model-dependent as it would assume both the SED of these sources 
(to predict their magnitudes at A <3 //m) and their z (to predict the location of their Lyman 
break and whether or not they are obser vable at A <3 /im). In any case, at 1.1 and 1.6 /im 



confusion is not reached until mAB~28 (iThompson et al 20051 ). If the first stars produced 



dusty environments their far-IR luminosities will be substantial and these sources should be 
visible at wavelengths redshifted today to mm and sub-mm bands. In that case, they may 
be resolvable with the ALMA0 large array, whose sub-mm resolution is better than 0.02". 

Finally, the fluctuations are unlikely to come from low-luminosity low- z normal galaxies. 
Such galaxies must have the surface density n > 3 x 10 7 deg -2 with the 3.6, 4.5 /im fluxes < 10-20 
nJy. Unless they are significantly fainter than this limit, emissions from star- forming systems 
should have comparable flux at shorter A out to the 4000 A break for passively evolving 
populations or to the Lyman cutoff at ~0.1 /im for star- forming gal axies. Galaxy counts now 



extend to m ~30.5 (2 nJy) at 0.67 /im and to 29 (10 nJy) at 1.6 /im (IMadau fc Pozzetti 20001 ) 
and are over an order of magnitude below the required value of n at the faintest magnitudes. 
This would exclude star forming galaxies as faint as 2 nJy at 0.67/(1 + z) /mi at z<5.7 and 
passively evolving populations out to 10 nJy at 1.6/(1 + z) /im at zt3. We note, however, 
that this analysis cannot exclude "abnormal" populations at low z. 

This work is supported by NSF AST-0406587 and NASA Spitzer NM0710076 grants. 
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Fig. 1. — Left: Shot noise power amplitude from the data is compared to the values of 
Psn estimated by i ntegrating the counts. Solid lines show the levels of Psn reached in the 
QSO1700 analysis flKAMMll ). Light shaded areas show the levels of P S n reached in KAMM2. 
Symbols plot Psn by integrating the counts evaluated for all five fields in Table 1 of KAMM2. 
Diamonds correspond to HDFN-E1 region, triangles to HDFN-E2, squares to CDFS-E1 and 
asterisks to CDFS-E2; open circles correspond to co unts for the QS 1700 field. Solid lin e 
shows Psn according to the fit to IRAC counts of (jFazio et al 2004T) used in (|KAMMlh ; 
dashed lines correspond to the IRAC counts analysis from (ISavage &: Oliver 20051 ). The 
counts are significantly incomplete due to confusion at the levels of Psn reached with our 
analysis and give a lower limit on the limiting magnitude. Right: CIB fluctuations from 
KAMM2 at the shot-noise levels shown with shaded regions in the left panel. The notations 
for the counts from the GOODS data is the same as in the left panels. Light-shaded dashes 
show the shot noise fluctuations. Solid lines show the least squares fits to the CIB fluctuations 
from sources at z>6 assuming ACDM model as described in the text. 
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Fig. 2. — Rest-frame luminosity per unit mass plotted vs wavelength for Population III 
spectra (from Santos et al - dashed lines) and "ordinary" stellar populations at 0.5 and 1 Gyr 
with Salpeter-Scalo IMF (computed from PEGASE for Z = 0, 10" 3 , 2 x 10" 3 , 5 x 10" 3 , KT 1 
assuming constant burst of star formation, i.e. SFR oc exp(— t/tburst) with tburst = 20 Gyr.) 
L & = 3.8 x 10 33 erg/sec is the solar bolometric luminosity. The part of emissions probed by 
the IRAC Channels 1 (3.6 jim) and 2 (4.5 jim) at z=5,10 is shown with the marked regions. 



